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ABSTRACT 

We present measurements of the anisotropic galaxy clustering from the Data Release 9 (DR9) 
CMASS sample of the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS). We an- 
alyze the broad-range shape of the monopole and quadrupole correlation functions to ob- 
tain constraints, at the effective redshift z=0.57 of the sample, on the Hubble expansion 
rate H(z), the angular-diameter distance Da(z), the normalized growth rate f(z)a%(z), the 
physical matter density il m h 2 , and the biased amplitude of matter fluctuation bas(z). We 
obtain {#(0.57), D A {0.57), /(0.57)er 8 (0.57), ft m h 2 , ba 8 (0.57)} = {87.6t^, 1396 ± 73, 

0.1261^-^0, 1-19 ± 0.14, 0.428 ± 0.066} and their covariance matrix as well. The param- 
eters which are not well constrained by our of galaxy clustering analysis are marginalized 
over with wide flat priors. Since no priors from other data sets (i.e., CMB) are adopted and 
no dark energy models are assumed, our results from BOSS CMASS galaxy clustering alone 
may be combined with other data sets, i.e. CMB, SNe, lensing or other galaxy clustering 
data to constrain the parameters of a given cosmological model. We show that the major 
power on constraining dark energy from the anisotropic galaxy clustering signal, as com- 
pared to the angular-averaged one (monopole), arises from including the normalized growth 
rate f(z)a%{z). In the case of the wCDM cosmological model our single-probe CMASS con- 
straints, combined with CMB (WMAP9+SPT), yield a value for the dark energy equation of 
state parameter of w — —0.90 ± 0.11. Therefore, it is important to include f(z)<j$(z) while 
investigating the nature of dark energy with current and upcoming large-scale galaxy surveys. 

Key words: cosmology: observations - distance scale - large-scale structure of Universe - 
cosmological parameters 
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1 INTRODUCTION 

The cosmic large-scale structure from galaxy redshift surveys 
provides a powerful probe of dark energy and the cosmologi- 
cal model that is highly complementary to the cosmic microwave 
background (CMB) (e,g.,[H inshaw et al. 2012), supernovae (SNe) 
dRiess et al.ll 19981: IPerlmutter et al.lfl999t) . and weak lensing (e.g. 
see lVan Waerbeke & MellierluOCOl for a review). 

The scope of galaxy redshift surveys has dramatically in- 
creased in the last decade. The 2dF Galaxy Redshift S urvey (2dF- 
GRS) obtai ned 221,414 galaxy redshifts at z < 0.3 dColless et al 



200 1L 12003). and the Sloan Digital Sky Survey (SPSS. I York et al. 
2000h has collected 930,000 galaxy spectra in t he Seventh Data 



Release (DR7) at z < 0.5 dAbazaiian etal]|2009l) . WiggleZ has 
collected spectra of 240,0 00 emission-line galaxie s at 0.5 < z < 1 
over 1000 square degrees (Drink water et al. Il2010t|parkinson et al.l 
2012b. and the Bary on Oscillation Sp ectroscopic Survey ( BOSS, 



Daws on etal.l2013r) of the SDSS-III (Eisen stein etalj201ll) is sur- 



veying 1.5 million luminous red galaxies (LRGs) at 0.1 < z < 0.7 
over 10,000 square degrees. The first BOSS data set has been made 
publi cly available recently in SDSS data release 9 (DR9. lAnn et al. I 
120121) . The planned space mission EuclidZ] will survey over 60 mil- 
lion emission-line galaxi es at 0.7 < z < 2 over 15,000 deg 2 
(e.g. lLaureiis et al.ll201lt) and the upcoming ground-based exper- 
iment Big BOS^] will survey 20 million galaxy redshifts up to 
z = 1.7 and 600,00 quasars (2.2 < z < 3.5) over 14,000 deg 2 
dSchlegeletal.l201ll) . The WFIRST satellite would map 17 million 
galaxies in the redshift range 1.3 < z < 2.7 over 3400 deg 2 , wit h 
a larger area possible with an extended mission dGreen et al.l2012h . 

Large-scale structure data from galaxy redshift surveys 
can be analyzed using either the power spectrum or the two- 
point correlation function. Although these two methods are 
Fourier transforms of one another, the analysis processes, the 
statistical uncertainties, and the systematics are quite differ- 
ent and the results cannot be converted using Fourier trans- 
form directly because of the finite size of the survey vol- 
ume. The SDSS-II LRG data have been analyzed, and the 
cosmologic al results delivered, using both t he power spectrum 
(see, e.g., |Tegmark etail 12004; Hutsi 2005; Padmanabha n et al.1 
2007 1; iBlake et alj |2007| ; iPercival et alj 120071. |2010| ; iReid et al.l 
20 id; iMontesano et al.1 120111). and the correlation function 



method (see, e.g., lEi senstein et lii] 20051; lOkumura et al.1 20081; 



Martinez etail 120091; 



'Sanchez et al.l 



Cabre & Gaztanagal 120091; 
20091; iKazinet alJ l2010d; Ichuang. Wang. & Hemanthal l2012i ; 
Samushia et al. 1 1201 lb IPadmanabhan et al.1 (20121 ; IXu et alj|2012h . 



Similar analysis have been also applied on the SDSS-III BOSS 
CMA SS sample and obtain ed the most p recise measurements to 
date dAnderson etal.112 012; Ma nera et aLll2013b iNuza et al.1120121 ; 
IReid et al.ll2012l ; ISamushia et al.ll2012l ; lToiei ro et al. 2Ql2). 

Galaxy clustering allows us to differentiate smooth dark en- 
ergy and modified gravity as the cause for cosmic acceleration 
through the simultaneous measurements of the cosmic expansion 
histor y H(z), and the g r owth rate of cosmic large scal e structure, 
f(z) Iguzzo et alJl2008l ; IWang|l2008l ; lBlake et alj|2012l) . However, 
to measure f(z), one must measure the galaxy bias b, which re- 
quir es measuring hig h er-order statistics of th e galaxy clustering 
(see lVerde et al.l2002l) . ISong & Percivail d2009l) proposed using the 
normalized growth rate, f(z )ag(z), which would avo id the uncer- 
tainties from the galaxy bias. lPercival & White! 12009) developed a 



1 http://sci.esa.int/euclid 

2 http://bigboss.lbl.gov/ 



metho d to measure f(z)ag(z) and applied it on simulations. IWand 
d2012l) estimated expected statistical constraints on dark energy 
and modified gravity, including redshift-space distortions and other 
constraints from galaxy clustering, using a Fisher matrix formal- 
ism. 

In principle, the Hubble expansion rate H(z), the angular- 
diameter distance Da{z), the normalized growth rate f(z)ag(z), 
and the physical matter density f2 m /i 2 can be well constrained by 
analyzing the galaxy clustering data alone. Eisenstein et al. ( 2005) 
demonstrated the feasibility of measuring Q m h 2 and an effective 
distance, Dv(z), from the SDSS DR3 LRGs, w here Dy(z) cor- 
respond s to a combination of H(z) and Da (z). Chuang & Wang 
(2012a) measured H(z) and Da{z) simultaneously using the 
galaxy clustering data from the two d imensional two-point corr e- 
lation function of SDSS DR7 LRGs. Ichuang & Wanal d2012bl lch 
improved the method and modeling to measure H(z), Da(z), /3, 
and ft m h 2 from the same d ata. 

ISamushia et al. I d201 II) mea sured f(z)ag(z) from the SDSS 
DR7 LRGs. iBlakeetalJ «2012h measured H(z), D A (z), and 
f(z)ag(z) from the WiggleZ Dark Energy Survey galaxy sample. 
IReid et al] J20 1 2) measured H( z), Da(z), and f(z)(jg, (z) from the 
SDSS BOSS DR9 CMASS and lSamushia et al) |2012|) derived the 
cosmological implications from these measurements to test devi- 
ations fr om the concordance A CDM model and general relativity 
(see also Nesse ris et al. I d20l 1) for using f(z)ag(z) to constrain 
modified gravity theories). 

In this study, we a pply the similar method and model as 
IChuang & Warigl d2012bl ld) to measure H(z), D A {z), f(z)a s (z), 
and Q, m h 2 which extracts a summary of the cosmological informa- 
tion from the large-scale structure of the SDSS BOSS DR9 CMASS 
alone. One can combine our single-probe measurements with other 
data sets (i.e. CMB, SNe, etc.) to constrain the cosmological pa- 
rameters of a given dark energy model. We also explore the strong 
power of adding f(z)ag(z) to the two dimensional galaxy cluster- 
ing analysis on constraining dark energy. 

This study is part of a series of papers performing anisotropic 
clustering analysis on the B OSS DR9 CMASS galaxy sample. 
Instead of using multipoles, Sanche z~et alj d2013h present a dif- 
fer ent method t aken f rom the 'clustering wedges' measurements 
by iKazin etal] J2013I) and combine the results with CMB, SNe, 
and Baryon Acoustic Oscillation s (BAO) to ob t ain co nstraints 
on the cosmological parameters. lAnderson et al.l (120 1 3h present 
anisotropic analysis using two approaches, multipoles and wedges, 
to obtain robust measurement of the BAO signal. 

This paper is organized as follows. In Section|2] we introduce 
the SDSS-III/BOSS DR9 galaxy sample and mock catalogues used 
in our study. In Section [3] we describe the details of the methodol- 
ogy that constrains cosmological parameters from our galaxy clus- 
tering analysis. In Section [4] we present our single-probe cosmo- 
logical measurements and demonstrate how to use our results with 
assuming different cosmological models or combining other data 
sets. In Section|5] we compare our results with previous or parallel 
works. In Section[6] we discuss the requirements to provide single- 
probe measurements. In Section|7] we apply some systematic tests 
to our measurements. We summarize and conclude in Sec. [8] 



2 DATA SET 

2.1 The CMASS Galaxy Sample 

The Sloan Dig i tal Sky Survey (SDSS: iFukugita et alj 1 19961: 



iGunn et alJl99l : lYork et al.l2000l;ISmee et al. 120121) mapped t 

© 0000 RAS, MNRAS 000, 000-000 



Single-Probe Cosmology Measurements from SDSS BOSS DR9 CM ASS 3 



one quarter of the sky using the dedicated 2.5m Sloan Telescope 
jAhn etal.1120121). The Baryon Oscillation Sky Survey (BOS S, 
lEisenstein et"aT]|201 1| ; iBolton et al. Il2012l : iDawson et al. |[2013h is 
part of the SDSS-III survey. It is collecting the spectra and red- 
shifts for 1.5 million galaxies, 160,000 quasars and 100,000 an- 
cillary targets. The Data Release 9 has been made public avail- 
ableJS We use galaxies from the SDSS-III BOSS DR9 CMASS cat- 
alogue in the redshift range 0.43 < z < 0.7. 'CMASS' samples 
are selected with an ap proximately constant stellar mass threshold 
jEisenstein et alj|201 lh . The sample we are using includes a total 
of 264,283 galaxies with 207,246 in the north and 57,037 in the 
south Galactic hemispheres. The median redshift of the sample is 
z = 0.57. The detai ls of generating this sample are described in 
lDawsonetaLU2013h . 



2.2 The Mock Catalogues 

Mane raet alj d2013l) created 600 mock catalogues for DR9 
CMASS sample. They created 2nd-order Lagrangian perturbation 
theory matter fields from which they populate haloes with mock 
galaxies using a halo occupation distribution prescription which 
has been calibrated to reproduce the clustering measure ments on 
scales between 30 and 80 ft _1 Mpc jwhite et al~ll201lh . We use 
these mock catalogues to construct the covariance matrix in our 
analysis. 



3 METHODOLOGY 

In this section, we describe the measurement of the multipoles of 
the correlation function from the observational data, construction 
of the theoretical prediction, and the likelihood analysis that leads 
to constraining cosmological parameters and dark energy. 



3.1 Measuring the Two-Dimensional Two-Point Correlation 
Function 

We convert the measured redshifts of the BOSS CMASS galaxies to 
comoving distances by assuming a fiducial model, i.e., flat ACDM 
with Q m = 0.274 and h = 0.7 w hich is the same mo del for con- 
structing the mock catalogues (see Manera et al. 2013) . We use the 
two-po int correlation function estimator given by lLandv & Szalavl 
dl993h : 



£(a,7r) = 



DD(a, tt) - 2DR(a, n) + RR(a, n) 



(1) 



RR(a, tt) 

where n is the separation along the light of sight (LOS) and a is 
the separation in the plane of the sky. DD, DR, and RR represent 
the normalized data-data, data-random, and random-random pair 
counts, respectively, for a given distance range. The LOS is de- 
fined as the direction from the observer to the center of a galaxy 
pair. The bin size we use is 1 /i _1 Mpc x 1 /i _1 Mpc. The Landy and 
Szalay estimator has minimal variance for a Poisson process. Ran- 
dom data are generated with the same radial and angular selection 
functions as the real data. One can reduce the shot noise due to ran- 
dom data by increasing the amount of random data. The number 
of random data we use is more than 15 times that of the real data. 
While calculating the pair counts, we assign to each data point a 
radial weight of 1/ [1 + n(z) ■ P w ] , where n(z) is the radial n umber 
density and P w = 2 ■ 10 4 h~ 3 Mpc 3 (see lFeldman et alJI 19941) . 

3 http://www.sdss3.org/ 



3.2 Theoretical Two-Dimensional Two-Point Correlation 
Function 

First, we adopt the cold dark matter model and the sim- 
plest inflation model (adiabatic initial condition). Thus, we can 
compute the linear matter power spectra, Pu n (k), by using 
CAMB (Code for Anisotropics i n the Microwave Background, 
I Lewis . Challinor, & Lasenbv 2000). The linear power spectrum can 
be decomposed into two parts: 



Plin{k)=Pun,(k)+P l BAo(k), 



(2) 



where P nw (k) is the "no-wiggl e" or pure CDM pow er spectrum 
calculated using Eq.(29) from lEisenstein & Hul d 19981) . P l BAo(k) 
is the wiggled part defined by the equation itself. The nonlinear 
damping effect of the "wig gled" part, in redshift space, c an be well 
approximated following lEisenstein. Seo. & White! ( 120071) by 



PBAC>{k,Hk) 



P l BAo{k)-exp 



2k': 



■[l + 4(2/ + / 2 )] ,(3) 



where /i^ is the cosine of the angle between k and the LOS, / is the 
growth rate, and fc* is co mputed following ICrocce & Scoccimarrol 
d2006l) ; lMatsubard d2008l) by 



1 

3^2 



Pun(k)dk 



-1/2 



The dewiggled power spectrum is 

Pd-w(k,^k) = Pn-w(k) + PSAo(k, /Xfe). 



(4) 



(5) 



Next, we include the linear redshift distortion as follows in 
order t o obtain the galaxy power spectrum in redshift space at large 
scales (Kaiser 1987|), i.e., 



p 9 fl (fc, Mfe ) = 6 2 (i + /?/4) 2 p d 4fc, Mfc ), 



(6) 



where b is the linear galaxy bias and f3 is the linear redshift distor- 
tion parameter. 

We compute the theoretical two-point correlation function, 
(,*{<y, 7r), by Fourier transforming the non-linear power spectrum 
Pg(k,fik)- This task is efficiently performed by using Legendre 
polynomial expa nsions and one-dimensi onal integral convolutions 
as introduced in Chuang&Wangld2012cl) . 

We convolve the 2D correlation function with the distribution 
function of random pair wise velocities , f v (v), to obtain the final 
model £ (a, it) following IPeeblesI dl980h by 



£(cr,7r) = / f <7,7r- 



H(z)a(z) 



fv(v)dv, 



(7) 



where the random mot ions (fingers of god) are represent ed by an 
exponential form (e.g jRatcliffe et al. Ill 9981 : lLandvll2002h 



,. , , 1 / y/2\v\ 

<7„V 2 



(8) 



where a v is the pairwise peculiar velocity dispersion. 

The cosmological parameter set that we use to compute the 
theoretical correlation function is {H(z), Da{z), /3(z), fl m h 2 , 
bas{z), Qbh , n s , cr v , f(z)}, where Q. m and fif, are the matter and 
baryons density fractions, n s is the power-law index of the primor- 
dial matter power spectrum, h is the dimensionless Hubble constant 
(Ho = WOh km s _1 Mpc _1 ), and cr&(z) is the normalization of the 
power spectrum. The linear redshift distortion parameter can be ex- 
pressed as P(z) = f(z)/b. Thus, one can derive f(z)as(z) from 
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the measured f3(z) and bas(z). On the scales we use for compari- 
son with the BOSS CMASS data, the theoretical correlation func- 
tion only depends on cosmic curvature and dark energy through the 
parameters H(z), Da(z), P(z), and b<7$(z) assuming that dark 
energy perturbations are unimportant (valid in the simplest dark 
energy models). Thus we are able to extract constraints from clus- 
tering data that are independent of a dark energy model and cosmic 
curvature. 



3.3 Effective Multipoles of the Correlation Function 

The traditional multipoles of the two-point correlation function, in 
redshift space, are defined by 

= ^ / dufcrfPtfa) 



where 



21 + 1 



V<r 2 + 7T 2 
COS /I, 



(9) 

(10) 
(11) 



and Pi (fi) is the Legendre Polynomial (I =0 and 2 here). We in- 
tegrate over a spherical shell with radius s, while actual measure- 
ments of £(<r, 7r) are done in discrete bins. To compare the mea- 
sured £(cr, tt) and our theoretical model, the last integral in Eq.(|9} 
should be converted into a sum. This leads to the definition for the 
effectiv e multipoles of the correlation function ( Chuang & Wang 
l2012bl) : 



E 



(2Z + l)C(<T,7r)P i (/i) v /l-^ 



Number of bins used in the numerator 
where As = 5 /i _1 Mpc in this work, and 

° = (n+ i)/i _1 Mpc,n = 0,1,2,... (13) 

tt = (m + i)/i _1 Mpc,m = 0, 1,2, ... (14) 

Both the measurement and the theoretical prediction for the ef- 
fective multipoles are computed using Eq.dl2K with f (a, 7r) given 
by the measured correlation function (see EqJT) for the measured 
effective multipoles, and Eq.lJTJ for the theoretical predictions. We 
do not use the conventional definitions of multipoles to extract pa- 
rameter constraints as they use continuous integrals (see Eq. [9^. 
Bias of the result could be introduced if the definitions of multi- 
poles differ between measurements from data and the theoretical 
model. 



3.4 Covariance Matrix 

We use the 600 mock catalogues created by iManera et al. 
for the BOSS CMASS DR9 to estimate the covariance matrix of 
the observed correlation function. We calculate the multipoles of 
the correlation functions of the mock catalogues and construct the 
covariance matrix as 



Cij — 



N 



(15) 



where TV is the number of the mock catalogues, X m is the mean of 
the m th element of the vector from the mock catalogue multipoles, 
and X^ is the value in the m th elements of the vector from the 
k th mock catalogue multipoles. The data vector X is defined by 
Eq.©. 

3.5 Likelihood 

The likelihood is taken to be proportional to exp(— x 2 /2) 
jPress et ail 19921) . with \ 2 given by 



X Qi 



bs,j\ 



(16) 



X = ^ ' [X t h,i — Xobsj] Cij 1 [X t h,j 

i J=l 

where Nx is the length of the vector used, Xth is the vector from 
the theoretical model, and X b s is the vector from the observed 
data. 

As explained in IChuang&Wan3 ( l2012j) . instead of recalcu- 
lating the observed correlation function while computing for differ- 
ent models, we rescale the theoretical correlation function to avoid 
rendering the \ 2 values arbitrary. It can be considered as an ap - 
plication of Alcock-Paczynski effect dAlcock & Paczvnskj|[l979l) . 
The rescaled theoretical correlation function is computed by 



where £th is computed by eq. I0, and \ 2 can b e rewrote as 

N x 

x 2 = {t 1 x t h,i — x^f^} Cfid,a ■ 



(17) 



{t- x x. 



th,j 



X 



obs,j 



} 



(18) 



where T~ X t h is the vector computed by eq. dl2t from the 
rescaled theoretical correlation function, eq. (1171 . X^ is the vec- 
tor from observed data measured with the fiducial model (see 
IChuang & W ang 2012a for more details regarding the rescaling 
method). 

3.6 Markov Chain Monte-Carlo Likelihood Analysis 

We perform Markov Chain Monte -Carlo likelihood analyses using 
CosmoMC jLewis & Bridldl2002h . The parameter space that we 
explore spans the parameter set of {_ff(0.57), Da(0.57), Q. m h 2 , 
(3(0.57), 6a 8 (0.57), Q b h 2 , n s , a v , /(0.57)}. Only {#(0.57), 
73.4(0.57), n m h 2 , /3(0.57), 6a 8 (0.57)} are well constrained us- 
ing the BOSS CMASS alone in the scale range of interest. We 
marginalize over the other parameters, {Qbh 2 , n s , <J V , /(0.57)}, 
with the flat priors {(0.01859, 0.02657), (0.865, 1.059), (0, 500) 
km/s, (0.5, 1)}, where the flat priors of Q, b h 2 and n s are cen- 
tered on the WMAP7 me asurements with a wi dth of ±7<twmap 
(aw map is taken from lKomatsu~et ai] |20ld) . These priors are 
sufficiently wide to ensure that CMB constraints are not dou- 
ble counted when our results ar e combined with CMB data 
dChuang, Wang. & Hemanth3l2012l) . 



4 RESULTS 

4.1 Measurement of multipoles 

Fig |l(a)| and |l(b)| show the effective monopole (£0) and quadrupole 
(£2) measured from the BOSS CMASS galaxy sample compared 
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with the theoretical model given the parameters measured. We are 
using the scale range, s — 40 — 160 /i _1 Mpc, and the bin size is 
5 /i _1 Mpc. The data points from the multipoles in the scale range 
considered are combined to form a vector, X , i.e., 



i So 



S2 ' 



t(2) 
?2 i 



;(»). 

! S2 i 



(19) 



where iV is the number of data points in each measured multipole; 
here N = 24. The length of the data vector X depends on the 
number of multipoles used. 



4.2 Measurement of Cosmological Parameters from BOSS 
CMASS only 

We now present the dark energy model independent measure- 
ments of the parameters {#(0.57), Da(0.57), tt m h 2 , 0(0.57), 
and bcr 8 (0.57)}, obtained by using the method described in 
previous sections. We also present the derived parameters in- 
cluding H (0.57) r s (z d ), D A {0.57)/r s (z d ), /(0.57)<r 8 (0.57), 
D v {0.57)/r s (z d ), A (0.57), a, and e 



D v (z) 



(l + z) 2 D A {zf 



H(z) 



A s {z) = D v (z) 



CZ 

Dv(z)/r s {z d ) 
{Dv[z)/r s {zd))fid ' 



and 



(H(z)D A (z)) ftd 



H(z)D A {z) 



1/3 



(20) 

(21) 
(22) 

(23) 



where r s (z d ) is the co moving sound hor i zon at the drag epoch cal- 
culated using eq. (6) in lEisenstein & Hul dl998h . 

Table Q] lists the mean, rms variance, and 68% confidence 
level limits for {#(0.57), L>a(0.57), tt m h 2 , /3(0.57), b<r 8 (0.57), 
#(0.57) r s (z d )/c, D A (0.57)/r s (z d ), D v (0.57)/r s (z d ), 
A s (0.57), a, and e} derived in an MCMC likelihood analysis from 
the measured £o + £2 of the DR9 CMASS correlation function. 

Table [2] gives the normalized covariance matrix for this pa- 
rameter set measured using £0 + £2- It is clear that the correla- 
tion between Q. m h 2 and Dv(0.57)/r s {z d ), D A (0.57)/r s (z d ), or 
r a (z d )H (0.57) /c are close to zero, since the dependency on Q, m h 2 
is removed by dividing or multiplying r s (z d ). 

For this measurement, we use 48 bins (£0 + £2), 9 fitting 
parameters (see Sec. 13. 6t . and scale range 40h~ Mpc < s < 
160ft _1 Mpc. The \ 2 P er degree of freedom (d.o.f.) is 0.51. This 
low value might indicate not only a good modeling but also possi- 
ble over-estimation for the covariance matrix constructed with the 
mock catalogues. 



Measured 


mean 


(j 


lower 


upper 


i7(0.57) 


87.6 


1.2 


80.8 


94.3 


-Da (0.57) 


1396 


74 


1324 


1470 


^ m h 


0.126 


0.019 


0.116 


0.134 


P = /(0.57)/6 


0.367 


0.084 


0.287 


0.446 


6o- 8 (0.57) 


1.19 


0.14 


1.05 


1.33 


Derived 


r s (z d )H(0.57)/c 


0.0454 


0.0031 


0.0426 


0.0482 


D A (0.57)/r s (z d ) 


8.95 


0.27 


8.69 


9.22 


D v (0.57)/r s (z d ) 


13.54 


0.29 


13.26 


13.82 


/(0.57)fr 8 (0.57) 


0.428 


0.069 


0.362 


0.494 


A s (0.57) 


0.436 


0.017 


0.419 


0.454 


a 


1.024 


0.022 


1.002 


1.045 


e 


0.015 


0.029 


-0.012 


0.042 



Table 1. The mean, standard deviation, and the 68% C.L. bounds 
of the measured parameters and the derived parameters from the 
BOSS DR9 CMASS galaxy clustering. The quantity a is defined as 
Dy(0.57) , D V (.0.57) fid _ ^ ^ Qf R ^ j^.j M -1 The ^ of 

r B (z d ) I r s (z d ) fid 

D A and r s (z d ) is Mpc. 



where 



_f/(0.57)- 8 7.6 \ 
J3a(0.57)-1396 ' 

n m h 2 -0.126 
6<t 8 (0.57)-1.19 
\ /(0.57)<T S (0.57)-0.42 8 / 



(25) 



and 



Ml 



\ 



0.03850 


-0.001141 


-13.53 


0.4007 


-1.271 


-0.001141 


0.0008662 


3.354 


-0.1598 


-0.3059 


-13.53 


3.354 


19370 


-987.8 


-770.0 


0.4007 


-0.1598 


-987.8 


110.8 


78.67 


-1.271 


-0.3059 


-770.0 


78.67 


411.3 



■ (26) 



/ 



One can use a subset of these parameters (measured 
and derived) and their covariance matrix to derive the cos- 
mological parameters. For example, if one is only interested 
in the cosmological parameters but not in the galaxy bias, 
b, one can use only four parameters, #(0.57), Da(0.57), 
Q, m h 2 , and /(0.57)<r 8 (0.57) to compute \ 2 t0 constrain 
the parameters of a given model. In Sec. 14.51 we use 
D v (0.57)/r s (z d ), {H(0.57)r s (z d ), LU(0.57)/r s (z d )}, and 
{H(0.57)r s (z d ), D A (0.57)/r 3 (z d ), /(0.57)cr 8 (0.57)} to explore 
the power on the constraining dark energy from /(0.57)<r 8 (0.57). 

In addition, we use H(z), D A (z), and Q m h 2 instead of 
H(z)r 3 (z d ), D A (z) /r s (z d ) to be more general. For example, 
while combining the supernovae data, which do not have Qi, h 2 as 
a parameter of the cosmological model, it is simpler to use H(z) 
than use H(z)r s (z d ). 

We also provide the code for using CosmoMC that includes 
BOSS CMASS clustering alonfl 



4.3 Using Our Results from CMASS only 

In this section, we describe the steps to combine our results with 
other data sets assuming some dark energy models. For a given 
model and cosmological parameters, including the linear galaxy 
bias b, one can compute #(0.57), D A (0.57), Q. m h 2 , bcr 8 (0.57), 
and /(0.57)u 8 (0.57). From Table Q] and [2] one can derive the co- 
variance matrix, Mjj, of these five parameters. Then, \ 2 can be 
computed by 

X 2 = A-cm^sM', 1 Acm A s.s, (24) 



4.4 Assuming Dark Energy Models 

In this section, we present examples of combining our CMASS- 
only clustering results with CMB data sets assuming specific dark 
energy models. 

Table [3] [4] [5] and[6]show the cosmological constraints assum- 
ing ACDM, oACDM (non-flat ACDM), wCDM (constant equation 

4 http ://members .iff. uam-csic . es/B OS SDR9singleprobe 
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Figure 1. Measurement of effective monopole (left) and quadrupole (right) of the correlation function for the BOSS DR9 CMASS galaxy sample (black 
points), compared to the the theoretical model given the parameters measured (solid line). The error bars are the square roots of the diagonal elements of the 
covariance matrix (see Sec. 13. 4) . In this study, our fitting scale range is 40/i — x Mpc < s < 160/i _1 Mpc. 





H 


D A 




P 


bug 


r a H/c 


if(0.57) 


1. 0000 


-0.2203 


0.5224 


0.1701 


0.1534 


0.8502 


£>a(0.57) 


-0.2203 


1.0000 


-0.7799 


0.4611 


-0.2938 


0.2530 




0.5224 


-0.7799 


1.0000 


-0.4100 


0.5746 


0.0237 


P = /(0.57)/6 


0.1701 


0.4611 


-0.4100 


1.0000 


-0.7220 


0.4219 


feo-g(0.57) 


0.1534 


-0.2938 


0.5746 


-0.7220 


1.0000 


-0.0952 


r s {z d ) ff(0.57)/c 


0.8502 


0.2530 


0.0237 


0.4219 


-0.0952 


1.0000 


D A (0.57) /r s {z d ) 


0.4038 


0.5739 


-0.0033 


0.3153 


0.1183 


0.4874 


D v (0.57) /r s (z d ) 


-0.5020 


0.2651 


-0.0200 


-0.1470 


0.2124 


-0.5828 


/(0.57)<j 8 (0.57) 


0.3565 


0.4387 


-0.1612 


0.8495 


-0.2768 


0.5300 


A s (0.57) 


0.0802 


-0.3945 


0.7669 


-0.4649 


0.6870 


-0.3090 


a 


-0.5020 


0.2651 


-0.0200 


-0.1470 


0.2124 


-0.5828 


£ 


-0.7989 


-0.4034 


-0.0102 


-0.4397 


0.0348 


-0.9476 




D A /r a 


D v /r s 




A a 


a 


e 


H(0.57) 


0.4038 


-0.5020 


0.3565 


0.0802 


-0.5020 


-0.7989 


£>a(0.57) 


0.5739 


0.2651 


0.4387 


-0.3945 


0.2651 


-0.4034 




-0.0033 


-0.0200 


-0.1612 


0.7669 


-0.0200 


-0.0102 


P = /(0.57)/6 


0.3153 


-0.1470 


0.8495 


-0.4649 


-0.1470 


-0.4397 


bo-g(0.57) 


0.1183 


0.2124 


-0.2768 


0.6870 


0.2124 


0.0348 


r s (z d )H(0.57)/c 


0.4874 


-0.5828 


0.5300 


-0.3090 


-0.5828 


-0.9476 


D A (0.57)/r a (z d ) 


1.0000 


0.4220 


0.5389 


0.2301 


0.4220 


-0.7316 


D v (0.57) /r s (z d ) 


0.4220 


1.0000 


-0.0515 


0.5431 


1.0000 


0.3081 


/(0.57)o- 8 (0.57) 


0.5389 


-0.0515 


1.0000 


-0.1512 


-0.0515 


-0.6026 


A s (0.57) 


0.2301 


0.5431 


-0.1512 


1.0000 


0.5431 


0.1678 


a 


0.4220 


1.0000 


-0.0515 


0.5431 


1.0000 


0.3081 


e 


-0.7316 


0.3081 


-0.6026 


0.1678 


0.3081 


1.0000 



Table 2. Normalized covariance matrix of the measured and derived parameters for the BOSS CMASS DR9 galaxy clustering. 



of state of dark energy), and otoCDM (non-flat universe with a con- 
stant equation of state of dark energy) models. In this study, we only 
list the parameters which can be well constrained by galaxy cluster- 
ing. We also present the results of the combination of CMASS and 
CMB data. The CMB data we use includes WMAP7 and WMAP9, 
which are the previous and the newest data release fr om the Wilkin- 
son Microwave Anisotropy Probe collaboration ; dKomatsu et al.l 
l20ld : iBennett et al J 12012 : Iffinshaw et al. Il2012h . We are also us- 
ing the newes t data release from the South Pole Telescope (SPT) 
collaboration dStorv et al. I [2012: iHou et al. Il20l3) . For WMAP7 
only and WMAP9 only data, we download the Markov chains from 



the WMAP websittQ' 0. While usin g WMAP9+SPT we o btain the 
Markov chains by using CosmoMC jLewis & B ridle 2002) with the 
data and likelihood c ode provide d by WMA P (Bennett et alJ2012r 
Iffinshaw et al. |2012» and SPT Istorv et al. |2012" : IHou et al. hod: 
Keisler et alj201 Th collaborations. 

One can see that the measurements from BOSS CMASS-only 
dataset are consistent with those from CMB, and adding CMASS 
to CMB produces significantly tighter constraints than using CMB 

5 WMAP7 http://lambda.gsfc.nasa.gov/product/map/dr4/parameters.cfm 

6 WMAP9 http://lambda.gsfc.nasa.gov/product/map/dr5/parameters.cfm 
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Figure 2. 2D marginalized contours for 68% and 95% confidence levels 
for O m and h (ACDM (left) or oACDM (right) model assumed) from 
WMAP9+SPT (red) and BOSS DR9 CMASS (green). The constraints from 
WMAP9+SPT are much better than those from CMASS only in ACDM 
model but are worse in the oACDM model. 



0.04 
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004 0.10 0.15 0.20 0.25 0.30 0.35 

LX m 

Figure 3. 2D marginalized contours for 68% and 95% confidence levels for 
n k and O m (oACDM model assumed) from WMAP9+SPT (red), CMASS 
(green), and WMAP9+SPT+CMASS (blue). The CMASS data break the 
degeneracy between and fi m constrained by CMB data. 



data alone. While adding SPT to WMAP9 in a ACDM model, fl m 
is decreased as found in IStorv et al. I (2012) (although they used 
WMAP7). It is interesting that the mean values from WMAP9+SPT 
in a ACDM model are much closer to those from WMAP7 than 
from WMAP9. 

Figure [2] shows that the constraints on f2 m and h obtained 
from CMB (WMAP9+SPT) are much better than CMASS-only in 
a ACDM model. However, the results are opposite in an oACDM 
model . Figure [3] shows how CMASS clustering breaks the degen- 
eracy between Q,k and Sl m constrained by CMB in the oACDM 
model, resulting in a much tighter constraint. Figure|4]demonstrates 
how CMASS clustering also breaks the degeneracy between w and 
Q. m constrained by CMB in the toCDM model, resulting in a much 
better constraint in which w is consistent (within 1 a) with w = — 1 
(cosmological constant model). This statement is true independent 
of which CMB data set (see Table [3J. Figure [5j shows that adding 
the CMASS and the CMB data improves the constraints on w and 
Qk significantly in the oiuCDM model, and the results are consis- 
tent with w = -1 and Q k = (i.e. a flat ACDM model). This 
statement holds regardless of which CMB data set is used (see Ta- 
ble [Hi. 



4.5 Power of the Constraints on dark energy from f(z)ag(z) 

In this section, we demonstrate how adding the measurement of 
f(z)cr&(z) from galaxy clustering data would improve the con - 
straints on the cosmological parameters. ISamushia et al. J 2012h. 
using the BOSS CMASS DR9 measurements from iReid et all 
d2012l) . found that the extra information from the 2D cor- 
relation function compared to the spherically-averaged corre- 
lation function improves the constraint on w significantly in 
the tuCDM model. The extra information from the anisotropic 
galaxy clustering includes the geometric distortion, also called 




0.40 



Figure 4. 2D marginalized contours for 68% and 95% confidence levels for 
w and fl m (u>CDM model assumed) from WMAP9+SPT (red), CMASS 
(green), and WMAP9+SPT+CMASS (blue). The CMASS data break the 
degeneracy between w and Q m constrained by CMB data. 



Alcock-Paczynski effect jAlcock & Paczvnskil 1 19791) . and the 
redshift space distortion (RSD). In Table. [7] we combine 
CMB (WMAP9+SPT) with different portions of the informa- 
tion obtained from our BOSS CMASS galaxy clustering anal- 
ysis. First, we use Dv{z)/r 3 , which is the main measure- 
ment from the spherically-averaged correlation function. Sec- 
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ACDM 




CMASS only 


WMAP7 only 


WMAP7+CMASS 


WMAP9 only 


WMAP9+CMASS 


WMAP9+SPT 


WMAP9+SPT+CM AS S 


Q m 


0. 


273 ± 0.032 


0.266 ± 0.029 


0.278 ± 0.019 


0.280 zb 0.026 


0.284 zb 0.017 


0.264 zb 0.019 


0.274 zb 0.015 


Hq 




68.0 ± 3.0 


71.0 zb 2.5 


69.8 zb 1.6 


70.0 zb 2.2 


69.5 zb 1.5 


71.2 zb 1.7 


70.3 zb 1.3 


6o- 8 (0.57) 




1.18 zb 0.14 


— 


1.18 ± 0.11 


— 


1.19 zb 0.10 


— 


1.17 zb 0.10 


/(0.57)cr 8 (0.57) 





.449 ± 0.055 


0.450 ± 0.025 


0.457 ± 0.018 


0.466 zb 0.019 


0.467 ± 0.015 


0.453 zb 0.016 


0.459 ± 0.012 


b 




2.01 zb 0.42 




1.94 zb 0.18 


- 


1.93 ± 0.17 




1.90 zb 0.17 


= /(0.57)/6 


0. 


387 ± 0.076 




0.391 ± 0.035 




0.396 zb 0.034 




0.397 ± 0.034 


n h 2 




zt U . Ulo 


u.ioo4t zt U.UUOO 


U.loDl zt U.UUoa 


U.-LoD^* zt U.UU'IO 


U.loOy zt U.UUoo 


U.looO zt U.UUOO 


U.lODZ zt U.UUoo 






0.80 zb 0.10 


0.800 zb 0.030 


0.806 zb 0.024 


0.821 ± 0.023 


0.820 zb 0.020 


0.806 zb 0.018 


0.811 ± 0.015 


/(0.57) 


0. 


750 ± 0.030 


0.743 zb 0.025 


0.753 ± 0.016 


0.754 =b 0.021 


0.758 zb 0.014 


0.742 zb 0.017 


0.750 ± 0.013 


if(0.57) 




90.3 ± 4.6 


94.2 ± 1.3 


93.52 =b 0.95 


93.9 zb 1.1 


93.59 zb 0.83 


94.36 zb 0.80 


93.96 ± 0.61 


D A (0.57) 




1413 ± 64 


1347 =b 33 


1364 ± 22 


1359 zb 29 


1366 zb 20 


1344 zb 21 


1355 ± 16 


D v (0.57) 




2106 zb 98 


2010 ± 42 


2031 ± 28 


2023 zb 37 


2033 zb 26 


2006 zb 27 


2019 ± 20 



Table 3. The cosmological constraints using different combinations of data assuming ACDM. One could see that the measurements from CMASS only are 
in good agreement with those from CMB. Combining CMASS with CMB gives significantly tighter constraints than using CMB only. The unit of H is 
kms -1 Mpc -1 . The unit of Da and Dy is Mpc. 



oACDM 


CMASS only 


WMAP7 only 


WMAP7+CMASS 


WMAP9 only 


WMAP9+CMASS 


WMAP9+SPT 


WMAP9+SPT+CMASS 




-0.05 ± 0.11 


-0.083 ± 0.082 


-0.0078 ± 0.0049 


-0.036 ± 0.059 


-0.0048 ± 0.0067 


0.005 ± 0.012 


-0.0043 ± 0.0054 




0.264 ± 0.033 


0.58 ± 0.30 


0.288 ± 0.019 


0.43 ± 0.23 


0.289 ± 0.022 


0.242 ± 0.052 


0.282 ± 0.019 


H 


68.8 ± 4.0 


53 ± 13 


67.5 ± 1.9 


62 ± 14 


68.5 ± 2.3 


75.7 ± 7.9 


68.9 ± 2.2 


bo- 8 (0.57) 


1.19 ± 0.14 




1.23 ± 0.11 




1.19 ± 0.12 




1.19 ± 0.11 


/(0.57)cr s (0.57) 


0.440 ± 0.058 


0.467 ± 0.028 


0.451 ± 0.020 


0.475 ± 0.029 


0.466 ± 0.018 


0.442 ± 0.024 


0.458 ± 0.012 


b 


2.10 ± 0.49 




2.09 ± 0.21 




1.96 ± 0.23 




1.98 ± 0.19 


= /(0.57)/6 


0.378 ± 0.076 




0.369 ± 0.036 




0.395 ± 0.042 




0.386 ± 0.035 


n m h 2 


0.125 ± 0.019 


0.1336 ± 0.0056 


0.1308 ± 0.0050 


0.1371 ± 0.0045 


0.1350 ± 0.0048 


0.1340 ± 0.0036 


0.1338 ± 0.0034 


°"8 


0.77 ± 0.11 


0.760 ± 0.046 


0.787 ± 0.030 


0.805 ± 0.036 


0.814 ± 0.026 


0.808 ± 0.020 


0.803 ± 0.018 


/(0.57) 


0.760 ± 0.040 


0.90 ± 0.12 


0.764 ± 0.016 


0.83 ± 0.11 


0.764 ± 0.018 


0.714 ± 0.053 


0.758 ± 0.016 


H(0.57) 


89.0 ± 5.0 


80.7 ± 9.8 


90.9 ± 1.8 


88 ± 10 


92.4 ±2.1 


98.3 ± 6.6 


92.5 ± 1.9 


Da (0.57) 


1408 ± 67 


1735 ± 309 


1407 ± 31 


1535 ± 268 


1385 ± 38 


1286 ± 110 


1379 ± 36 


£>v(0.57) 


2112 ± 101 


2516 ± 397 


2094 ± 44 


2252 ± 349 


2061 ± 53 


1924 ± 153 


2054 ± 50 


Table 4. Same as 


Table[3]but assuming 


; oACDM. We note that the constraints on the curvature, Qj., 


is in a good agreement with a flat universe (O^ = 0). 


iuCDM 


CMASS only 


WMAP7 only 


WMAP7+CMASS 


WMAP9 only 


WMAP9+CMASS 


WMAP9+SPT 


WMAP9+SPT+CMASS 


w 


-1.24 ± 0.42 


-1.09 ± 0.38 


-0.94 ± 0.12 


-1.01 ± 0.43 


-0.94 ± 0.13 


-0.97 ± 0.37 


-0.90 ± 0.11 


n m 


0.247 ± 0.064 


0.264 ± 0.098 


0.290 ± 0.029 


0.30 ± 0.11 


0.295 ± 0.034 


0.30 ± 0.11 


0.297 ± 0.031 


H 


73.0 ± 9.2 


75 ± 13 


68.3 ± 3.4 


71 ± 14 


68.1 ± 3.8 


71 ± 13 


67.5 ± 3.5 


6<r 8 (0.57) 


1.20 ± 0.14 




1.20 ± 0.11 




1.19 ± 0.10 




1.21 ± 0.11 


/(0.57) CTS (0.57) 


0.435 ± 0.065 


0.471 ± 0.091 


0.442 ± 0.033 


0.47 ± 0.10 


0.450 ± 0.035 


0.450 ± 0.079 


0.436 ± 0.029 


b 


2.20 ± 0.58 




2.04 ± 0.26 




2.00 ± 0.25 




2.08 ± 0.26 


= /(0.57)/b 


0.370 ± 0.082 




0.372 ± 0.049 




0.381 ± 0.050 




0.364 ± 0.048 


n m h 2 


0.127 ± 0.019 


0.1335 ± 0.0056 


0.1342 ± 0.0039 


0.1364 ± 0.0048 


0.1357 ± 0.0037 


0.1339 ± 0.0036 


0.1343 ± 0.0030 


"a 


0.75 ± 0.15 


0.82 ± 0.13 


0.785 ± 0.046 


0.82 ± 0.15 


0.796 ± 0.048 


0.79 ± 0.12 


0.778 ± 0.041 


/(0.57) 


0.773 ± 0.055 


0.752 ± 0.029 


0.747 ± 0.016 


0.764 ± 0.027 


0.752 ± 0.017 


0.750 ± 0.020 


0.745 ± 0.014 


ff(0.57) 


88.4 ± 6.1 


93.1 ± 1.9 


93.8 ± 1.0 


92.6 ± 1.9 


93.72 ± 0.91 


93.3 ± 1.4 


93.99 ± 0.61 


D A (0.57) 


1400 ± 69 


1346 ± 104 


1375 ± 30 


1384 ± 115 


1377 ± 33 


1374 ± 110 


1379 ± 31 


£V(0.57) 


2109 ± 101 


2015 ± 108 


2040 ± 33 


2057 ± 119 


2043 ± 35 


2042 ± 114 


2043 ± 33 


Table 5. Same as Table[3]but assuming 


; tuCDM. One can see that the constraints 


on w is consistent with -1 (cosmological constant model). 




oiuCDM 


CMASS only 


WMAP7 only 


WMAP7+CMASS 


WMAP9 only 


WMAP9+CMASS 


WMAP9+SPT 


WMAP9+SPT+CMASS 


w 


-1.19 ± 0.45 


-1.18 ± 0.60 


-0.99 ± 0.21 


-0.88 ± 0.59 


-0.93 ± 0.14 


-0.77 ± 0.35 


-0.91 ± 0.13 




-0.02 ± 0.13 


-0.105 ± 0.084 


-0.0017 ± 0.0084 


-0.053 ± 0.063 


-0.0010 ± 0.0067 


-0.009 ± 0.025 


-0.0014 ± 0.0071 


n m 


0.258 ± 0.061 


0.67 ± 0.32 


0.297 ± 0.044 


0.56 ± 0.25 


0.302 ± 0.033 


0.38 ± 0.18 


0.299 ± 0.030 


H 


70.9 ± 8.3 


50 ± 16 


67.9 ± 4.8 


54 ± 14 


67.4 ± 3.7 


65 ± 16 


67.1 ± 3.3 


6<t 8 (0.57) 


1.19 ± 0.14 




1.27 ± 0.11 




1.21 ± 0.11 




1.22 ± 0.11 


/(0.57) CT8 (0.57) 


0.442 ± 0.063 


0.484 ± 0.077 


0.454 ± 0.042 


0.45 ± 0.10 


0.450 ± 0.039 


0.412 ± 0.066 


0.436 ± 0.032 


b 


2.12 ± 0.57 




2.16 ± 0.28 




2.04 ± 0.26 




2.11 ± 0.25 


= /(0.57)/6 


0.380 ± 0.081 




0.361 ± 0.058 




0.376 ± 0.049 




0.360 ± 0.047 


n m h 2 


0.126 ± 0.019 


0.1345 ± 0.0056 


0.1346 ± 0.0048 


0.1365 ± 0.0048 


0.1360 ± 0.0039 


0.1339 ± 0.0037 


0.1338 ± 0.0035 




0.78 ± 0.14 


0.76 ± 0.11 


0.793 ± 0.051 


0.74 ± 0.14 


0.794 ± 0.050 


0.73 ± 0.11 


0.775 ± 0.042 


/(0.57) 


0.762 ± 0.056 


0.94 ± 0.11 


0.764 ± 0.027 


0.87 ± 0.10 


0.756 ± 0.025 


0.760 ± 0.064 


0.748 ± 0.021 


H(0.57) 


89.2 ± 6.0 




92.1 ± 3.6 




93.4 ± 2.5 




93.5 ± 2.4 


D A (0.57) 


1406 ± 69 




1396 ± 41 




1387 ± 41 




1388 ± 35 


£>v(0.57) 


2108 ± 101 




2075 ± 58 




2055 ± 54 




2056 ± 47 



Table 6. Same as Table[3]but assuming omjCDM. In all the cases, the constraints are consistent with w = — 1 and Slfc = (ACDM model). 
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-0.04 -0.02 0.00 0.02 0.04 

Figure 5. 2D marginalized contours for 68% and 95% confidence levels 
for w and Q fc (oiuCDM model assumed) from WMAP9+SPT+CMASS 
(blue). While w and are not well constrained by WMAP9+SPT or 
CMASS-only, we only show the contour of 68% confidence level from 
WMAP9+SPT (red). Adding the CMASS data on the CMB data improves 
the constraints on w and significantly, and the results are consistent with 
w = -1 and tt k = (ACDM model). 

ond, we use H(z)r 3 +DA(z)/r 3 , which could be considered as 
adding geometric distortion on the previous one. Third, we use 
H(z)r a +DA(z)/r s +f(z)ag(z), which consists of adding the RSD 
information. We find that adding f(z)a%{z) actually dominates the 
improvement of the cosmological parameter constraints; but there 
is only a small difference between using H(z)r 3 +DA(z)/r„ and 
using Dv{z)/r s . There is no significant improvement between 
using H(z)r s +DA( z)/r s and using Dy(z)/r s measured from 
SPSS DR7 LRGs JWang, Chuang & Mukherieel 1201 j ; IXu et alj 
|2012|) . lAnderson et alj d2013l) also find the similar results from 
BOSS DR9 CMASS galaxy sample. 

Figure [6] shows the constraints obtained from using different 
portions of the information from galaxy clustering as described 
above. There is not much difference when replacing Dv(z)/r„ 
with H(z)r s and Da{z)/t s , but significant improvem ent when 
adding f(z)ag(z). As explained in Samush ia et alj J2012h . the cor- 
relation between f(z)as(z) and the geometric distortion (Alcock- 
Paczynski effect) contributes to the improvement of the constraint 
on dark energy as well. 



5 COMPARISON WITH PREVIOUS WORKS 

We compare our results with previous or parallel works which use 
the same data (BOSS CMASS DR9). In general, although the data 
set used are identical, the results could be slightly different due 
to using different methodology, models, or parameter spaces ex- 
plorer! 

lAnderson etail d2012l) measured a = 1.016 ± 0.017 from 
the position of the BAO peak of the two-point correlation function 
(monopo le), before applying th e reconstruction method on the den- 
sity field. lSanchez et alj d2012h measured a = 1.015 ± 0.019 from 



CMB + D v lr, 
CMB + Hr, + D A /r s 
CMB + Hr, + D A /r s + f-cr s 




0.35 0.40 0.45 0.50 0.55 0.60 0.65 
f(0.57)cr 8 (0.57) 



Figure 6. 2D marginalized contours for 68% and 95% confidence levels 
for w and /(z)crg(z) (uiCDM model assumed) from CMB+Dv(z)/r B 
(thin solid red), CMB+.H"(z)r s +D A (z)/r' s CMASS(dashed blue), and 
CMB+H(z)r 3 +DA(z)/r s +f(z)as(z) (thick solid black). One can see 
that there is no much improvement while replacing Dy(z)/r a with 
H(z)r s +DA(z)/r s , but have significantly improvement on dark energy 
while adding f(z)crg(z). 



the fu ll shape of the correlation function (monopole). iReid et all 
(2012) obtain a — 1.023 ± 0.019 from the broad-range shape 
of the monopole and quadmpole of the correlation function; and 
iRoss et al.1 ( f2012h obtain a = 1.020 ± 0.019 using monopole and 
quadrupole as well. These results are all consistent with our result 
a = 1.024 ± 0.022. Our uncertainty is slightly larger because we 
marginalize over wide ranges of Qbh 2 and n s as mentioned in Sec 

m 

Comparing to our parallel works, lAnderson et al.l d20 1 3b ob- 
tain £> A (0.57)/r s = 9.20 ± 0.29 and #(0.57)r s /c = 0.0474 ± 
0.0040; iKazin et alj J2013t) obtain Da (0.57) /r s = 9.05 ± 0.27 
and H(0.57)r s /c = 0.0464 ± 0.0030; both are in excellent agree- 
ment with our measurements, Z?A(0.57)/r s = 8.95 ± 0.27 and 
H(0.57)r a /c = 0.0454±0.0031, despite of the differences amon g 
different methodology, models, and analysis. ISanchez et al.N2013h . 
our other parallel work, combines CMASS and other datasets 
(CMB, SNe, etc) to measure cosmological parameters adopting 
specific dark energy models. Here, we only discuss the cases 
of combining CMASS and CMB datasets for comparison. When 
combining CMASS and CMB for ACDM, they obtain fi m = 
0.285 ± 0.015 and Ho = 69.5 ± 1.2, which agrees to our results 
Q m = 0.274 ± 0.015 and H = 70.3 ± 1.3. In the case of a more 
general model, i.e. okiCDM, they obtain fi fc = -0.0023 ± 0.0061 
and w — —0.97 ± 0.16, which results are consistent with ours 
n fe = -0.0014 ± 0.0071 and w = -0.91 ± 0.13. Note, that de- 
spite of different model, analysis, and CMB dataset, their results 
are i n good agre e ment with ours. 

IReid et alj J2012h obtained {/(0.57)a 8 (0.57), F(0.57), 
D v (zQ.57)/r s } ={0.4298 ± 0.0672, 0.6771 ± 0.0417, 1.0227 ± 
0.0188} where F(z) = (1 + z)D A (z)H(z)/c. Our corresponding 
measurements are {0.428±0.069, 0.637±0.057, 1.024±0.022}. 
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0.1340 ± 0.0036 


0.1342 ± 0.0036 


^8 


0.751 ± 0.063 


0.772 =b 0.060 


0.780 ± 0.041 


0.752 =b 0.084 


0.779 ± 0.081 


0.782 ib 0.045 


/(0.57) 


0.746 ± 0.016 


0.747 =b 0.017 


0.748 ± 0.014 


0.749 ib 0.029 


0.753 ± 0.031 


0.753 ib 0.021 


Jf(0.57) 


93.68 ± 0.70 


93.81 ± 0.70 


93.91 ± 0.61 


93.3 ib 2.8 


93.0 ib 3.0 


93.1 ib 2.3 


D A (0.57) 


1407 ± 44 


1389 ± 37 


1382 ± 28 


1410 ib 43 


1392 ± 36 


1390 ± 30 


D v (0.57) 


2073 ± 45 


2054 =b 38 


2047 ± 29 


2078 ib 46 


2064 ± 42 


2060 ± 40 



Table 7. The cosmological constraints using CMB (WMAP9+SPT) and different portions of the information extracted from CMASS assuming tuCDM and 
otoCDM. The unit of\ff is kms - 1 Mpc -1 and the unit of Da and Dy is Mpc. 



F(z), in our analysis, can be derived from e, see Eq. l |23l l. Although 
the results are in good agreement, their 6.2% error on F(z) is sig- 
nificantly smaller than ours (8.9%). The reason should be that they 
included smaller scales in their analysis, thus they could obtain 
tighter constraints. While the uncertainties from galaxy bias, red- 
shift space distortion, and nonlinear effect have larger effects on 
small scales, we use the scale range which is not sensitive to these 
uncertainties to obtain robust cosmological measurements from the 

galaxy sample. 

Regarding other results, iBlake et aL I (120121) measured 
{/(z)o-g(z), H(z), Da(z)} = {0.390 ± 0.063, 87.9 ± 6.1, 
1380 ± 95} at z = 0.6 usi ng the galaxy sample from Wig- 
gleZ dDrinkwater et al. I l2010h . which is consistent with our 
measurements {0.428 ± 0.069, 87.6 ± 7.2, 1396 ± 74} at 
z = 0.57. 



6 SINGLE-PROBE MEASUREMENTS: REQUIREMENTS 

In general, an ideal single-probe cosmological measurement should 
not adopt any cosmological parameter priors from another data set 
(e.g., CMB, SNe, ...). Otherwise, one should not be able to combine 
the results with other cosmological probe data set to avoid double 
counting. However, in practice, some priors are needed because not 
all the cosmological parameters could be well constrained by the 
given observed data. In other words, the data set can not constrain 
all the parameters of a given model, so that we need to use priors 
on some of them. In addition, adopting priors can reduce computing 
time significantly, which might be important when using a compli- 
cate theoretical model. In this section, we discuss how priors may 
be chosen in order to avoid biases in the measurements and under- 
estimations of the uncertainties. 

Which cosmological parameters are appropriated to be set on 
the priors? The first choice would be the parameters which are 
tightly constrained by other data sets. For example, the uncertain- 
ties of n s and Qbh 2 measured by CMB could be lower than 1.5% 
and 2.5% respectively. Therefore, the priors on these two parame- 
ters are often used in the analyses of galaxy clustering data sets. 

What prior we should use for a given parameter? The key is 
that we have to make sure that the prior will not result in double 
counting when later we combine with other probes. This can be 
done by using a sufficiently wide flat prior. For example, in this 
study we use ±7<jwmap for n s and Qbh 2 ; where owmap is the 
uncertainty measured from WMAP7. 

Due to some practical concerns (e.g., the data set size or com- 
puting time), one might want to use tighter priors or to adopting 
priors for more parameters. It is fine if the quantities one is mea- 



suring are insensitive to those parameters with priors. For example, 
as shown in Table [2] the correlation between Dv(z)/r s (Zd) and 
f2 m /i 2 is close to zero (i.e. -0.0200), so that one could put a strong 
prior on Q. m h 2 (i.e., CMB prior) and measure Dy(z) /r a (zd) with- 
out introducing a bias or underestimation. Yet, the proper quanti- 
ties to be mea sured might vary w ith different analysis conditions. 
For example, iBlack et al] d201 lal) found that A s (z), instead of 
Dv(z)/r s (zd), is uncorrected to fi m /i 2 while including smaller 
scales (see Fig. 6 and 15 in their paper). In other words, when 
including small scales, Dv(z)/r s (zd) is no longer a good mea- 
surement while u sing Sl m /i 2 prior. It is worth to mention that 
lBlacketalJd2011iJ) marginalized over Q m h 2 , which is similar as 
we did in our study, so that both A s (z) and Dv{z) /r 3 (zd) are fine 
measurements in their study (they recommended to use A s (z) as 
their default result). 

In this study, we provide the single-probe clustering measure- 
ments by using sufficient wide priors on the parameters, which are 
not well constrained by BOSS CMASS data. In addition, it is worth 
to emphasize that our results are obtained without assuming a dark 
energy model. Thus, one can derive cosmological parameter con- 
straints of a specific dark energy model from our results or combin- 
ing our results with other datasets (see examples in Sec. 14.41 . 



7 SYSTEMATIC TESTS 

iRoss et al.1 j2012h studied the systematics of the BOSS CMASS 
data selection, observation, selection functions of the geometry of 
the survey. In this section, we focus on the possible systematics 
from the theoretical model due to the scale-dependent uncertain- 
ties including nonlinear effects, galaxy bias and redshift distortions. 
These effects mainly affect the small scales and they might actually 
be cancelled by each oth er. For example, t he galaxy bias tends to be 
larger at smaller scales dNuza et al.ll2012l) , but the nonlinear effect 
makes t he correlation function smaller at t he same scale (e.g. see 
Fig 3 in IChuang, Wang. & Hemanthll2012l) . Therefore, if one only 
corrects the model with one of the effects (i.e. nonlinear effect), 
the corrected model could become less accurate. Here, instead of 
adding corrections, we use the scale range which is not sensitive 
to these uncertainties to obtain robust cosmological measurements 
from the galaxy sample. Table [8] shows the systematic tests that 
we have performed by varying the scale range used. We vary the 
lower limit from 25 to 50 ft _1 Mpc and the upper limit from 120 
to 200 /i _1 Mpc. We find that the results are insensitive to those 
scale limits, which demonstrates that our constraints are robust to 
the combination of scale-dependent uncertainties mentioned above. 
We also apply our method on the BOSS CMASS sample in the 
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north Galactic hemispheres and south Galactic hemispheres sepa- 
rately. We find that a = 0.996^ g|| and 1.085± ^g|g for north and 
south, which is in excellent agreement with what lRoss et 
obtained, i.e., a = 0.994 ± 0.023 and 1.083 ± 0.029. While the 
deviat ions between the results from north and south samples are 
large, iRoss et al.l (2012) found that it is probably due to the statisti- 
cal variance instead of the presence of some systematic error. 



8 SUMMARY 

In this paper, we present measurements of the SDSS-III/BOSS 
DR9 anisotropic galaxy clustering and provide cosmological 
constraints from our CMASS data set only, and in combination 
with other investigations in particular with CMB (WMAP9 and 
SPT) data. We summarize our study as follows: 

(i) We present single-probe measurements of H(0.57), 
D A (0.57), /(0.57)a 8 (0.57), and Q m h 2 as a summary of the 
information extracted from the BOSS CMASS galaxy clustering 
signal (results are listed in Table Q] and |2}. A CosmoMC code that 
includes our BOSS CMASS clustering alone is providecQ. 

(ii) Our cosmological constraints are obtained without 
assuming a flat universe or dark energy model. One can combine 
our results with other probe data sets and derive cosmological 
implications for a given model. We have explained the steps 
required to use our results and also demonstrated its applications 
with some examples. 

(iii) Combining our results from CMASS-only with CMB 
data, we find that the constraint on the constant equation of 
state of dark energy, w, can be significant improved by adding 
the measurement of f(z)ag(z) to the H(z)r s and DA(z)/r 3 
measurements, i.e., f(z)ag(z) provides the strong power of 
the improvement going from spherical averaged analysis to the 
anisotropic galaxy clustering analysis. Of course, the correlation 
among all the measu rements must be taken into account; see 
ISamushia etail d2012h . Our results are all consistent with w = — 1, 
which corresponds to the cosmological constant model (i.e. 
ACDM). 

(iv) We have compared our results with other investigations 
using the same BOSS DR9 CMASS data and find all the previous 
and new results are in excellent agreement. We also found that our 
results are insensitive to the scale range used, which is likely due to 
the cancellation of the nonlinear effect and scale-dependent bias. 

(v) We explain and discuss the requirements of a given anal- 
ysis to provide a single-probe measurement (see Sec. [6} in order 
to avoid double counting when combine with other cosmological 
probe data sets. 

Our methodology can be applied on the current and future 
large-scale galaxy surveys (e.g. eBOSS, BigBOSS, and Euclid) 
to obtain single-probe and model independent cosmological con- 
straints, which will provide a powerful and convenient way to per- 
form a joint data analysis with other data sets. 
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Table 8. This table presents the systematic tests with the scale range and the regime of the sample used. The fiducial results are obtained by considering the 
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